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A High Statistics Lattice Calculation of Quark Masses with a 
Non-Perturbative Renormalization Procedure. 

L. GiustP* 

'^Scuola Normale Superiore and INFN Sczionc di Pisa, P. zza dei Cavalieri 7, 1-56100 Pisa Italy 

We present results of a high statistics study (0(2000) configurations) of the quark masses in the MS scheme 
from Lattice QCD in the quenched approximation at /3 = 6.0, fH — 6.2 and /3 = 6.4 using both the Wilson and the 
tree-level improved SW-Clover fermion action. We extract quark masses from the meson spectroscopy and from 
the axial Ward Identity using non-perturbative values of the renormalization constants. We compare the results 
obtained with the two methods and we study the 0{a) dependence of the quark masses for both actions. Our 

best results are mf ^(2 GeV) = (123 ± 4 ± 15) MeV and ^(2 GeV) = (1525 ± 40 ± 100) MeV. 



1. INTRODUCTION 

Quark masses are fundamental parameters of 
the QCD Lagrangian that cannot be measured 
directly in the experiments, since free quarks are 
not physical states, and that cannot be predicted 
by QCD using only theoretical considerations. 
Up to now QCD sum rules and lattice simulations 
are the only non-perturbative techniques able to 
determine the absolute values of the light quark 
masses. Lattice technique does not require any 
additional model parameter and is the only pro- 
cedure that can be systematically improved. 

2. QUARK MASSES 

The usual on-shell mass definition cannot be 
used for quark masses and their values depend 
on the theoretical definition adopted. In the fol- 
lowing we will give our final results for the quark 



masses defined in the AIS scheme. 



(^) is de- 



fined by the perturbative expansion of the quark 
propagator renormalized with the MS prescrip- 
tion [Q and it depends on the renormalization 
scale /X. The symbols used in the following are 
fully defined in |l|,|. 

The "bare" lattice quark mass can be determined 
directly from lattice simulations by fixing the 
mass of a hadron containing a quark with the 
same flavour to its experimental value. 
Quark masses can also be extracted from the ax- 
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ial Ward Identity. Close to the chiral limit and 
neglecting terms of 0{a), the Ward Identity can 
be written as 
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where Za, Zp and Zs are the renormalization 
constants of the axial, pseudoscalar and scalar 
densities and nic is defined in |l|,|. The stan- 
dard perturbative approach uses the lattice and 
the continuum perturbation theory to connect 
the "bare" lattice quark mass to the m^^^{fi) 
The scale 1/a, where a is the lattice spac- 
ing, of our simulations is a^^ ~ 2 — 4 GeV. At 
these scales we expect small non-perturbative ef- 
fects. However the "tadpole" diagrams, which are 
present in the lattice perturbation theory, can 
give raise to large perturbative corrections and 
then to large uncertainties in the matching pro- 
cedure at values of /3 = = 6.0 — 6.4. 
The non-perturbative renormalization (NP) tech- 
niques eliminate these uncertainties [0lq|. Za, Zp 
and Zs can be calculated by imposing the renor- 
malization conditions, described in the next sec- 
tion, on the quark states of momentum — fi^ 
and in the Landau gauged]. The quark mass in 
the M S scheme is then defined as 
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Table 1 



Summary of the parameters 


of the runs analyzed in this work. 














Matrix Elements 








C60a 


C60b 


C60c 


C60d 


W60 


C62a 


W62a 


W62b C64 


W64 


P 6.0 


6.0 


6.0 


6.0 


6.0 


6.2 


6.2 


6.2 6.4 


6.4 


Action SW 


SW 


SW 


SW 


Wil 


SW 


Wil 


Wil SW 


Wil 


# Confs 490 


600 


200 


200 


320 


250 


250 


110 400 


400 


Volume 18^ x 64 24^ x 40 18^ x 32 


16^ X 32 


18^ X 


64 243 X 64 243 X 


64 243 X 64 243 X 


64 243 X 64 










Renormalization Constants 












C60Z 


W60Z 


C62Z 


W62Z 


C64Z 


W64Z 








6.0 


6.0 


6.2 


6.2 


6.4 


6.4 




Action 




SW 


Wil 


SW 


Wil 


SW 


Wil 




# Confs 




100 


100 


180 


100 


60 


60 




Volume 




16^ X 32 


16^ X 32 


163 X 32 163 X 32 


243 X 


32 243 X 32 





Table 2 

Quark Masses from the spectroscopy in MeV. MS 



masses are at ji 


= 2 GeV 






Run 






TOc(a) 




C60a 


83(2) 


118(7) 






C60b 


81(2) 


117(7) 






C60c 


83(3) 


119(7) 






C60d 


79(3) 


113(7) 






W60 


98(2) 


133(16) 


1335(31) 


1612(316) 


C62a 


83(4) 


120(9) 


1106(49) 


1470(170) 


W62a 


93(3) 


131(15) 


1205(24) 


1553(249) 


W62b 


92(4) 


129(16) 


1206(32) 


1560(257) 


W64 


82(3) 


120(16) 






C64 


69(3) 


103(9) 







where 

C/r^(M,M') = ( 



MS(,, ,,r\ _ ( O'^(m) ■|7'°'/2/3o 



a,(M)-a,(A.O 7W 7(°)/3i. 
47r ^ 2/3o 



and, for large time separations, 

1 (Ao(r)P(O)) 
^("^ = 2^'^"' {Pir)Pm ■ 



(3) 



(4) 



This procedure uses only the continuum pertur- 
bation theory to connect the quark mass from 
the i?/-scheme to MS. The continuum pertur- 
bation theory is used at scales /i ~ 2 — 4 GeV 



large enough to avoid non-perturbative effects or 
higher order corrections. 

The non-perturbative approach can be applied di- 
rectly to the unquenched case 

3. NP RENORMALIZATION 

For a generic two-quark operator Or (a) = 
iliTip^ we define the renormalized operator 
Ot{ij) = ZQ^Or{a) by introducing the renormal- 
ized constant Zq^ calculated imposing the renor- 
malization condition [Q 



zRI/ 



1 . 



^Or(Ma)^^ ifia)Tor{pa)\p2=f,2 

(po-) is the forward amputated green func- 
tion of the bare operator Or (a) on off-shell quark 
states with = fi'^ in the Landau gauge. For 
we have used a definition inspired by the Ward 
Identities 



Z,h — 



Tr (j2x=iA 7A sin{pxa)S ^ (pa)) 



where S{pa) is the quark propagator. These 
renormalization conditions satisfy the continuum 
Ward Identity, i.e. Zp/Zs = Zj^^/Z§^. 
This procedure works if fi satisfies the condi- 
tion Aqcd /X <C 1/a to avoid large higher- 
order perturbative corrections and discretization 
errors. In the figure |l| we report the matrix el- 
ement of the renormalized pseudoscalar operator 
in the chiral limit as a function of the renormal- 
ization scale for the runs we have analysed. Data 
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Table 3 

Quark Masses from the Ward Identity in MeV. 
M S masses are at /i = 2 GeV. 



Run 


Ps{a) 




Pc{a) 




C60a 


57(2) 


129(6) 






C60b 


55(1) 


125(5) 






C60c 


59(2) 


132(6) 






C60d 


55(2) 


127(6) 






W60 


77(2) 


126(4) 


956(20) 


1557(52) 


C62a 


60(4) 


127(9) 


747(25) 


1572(71) 


W62a 


75(3) 


120(6) 


936(16) 


1502(47) 


W62b 


71(3) 


115(6) 


927(22) 


1491(55) 


W64 


69(4) 


107(7) 






C64 


55(4) 


106(8) 








Figure 1. Renormalte^tfySSkdoscalar density ver- 
sus log(/x^) for the runs analysed. 



show that the discretization errors are within the 
statistical errors in the range /io ~ 1 where the 
NP-renormalization is applied to rcnormalize the 
quark masses. All details of NP-rcnormalization 
will be given in a forthcoming paper ||]. 

4. RESULTS 

We have calculated the quark masses from the 
meson spectroscopy and from the Ward Identity 
using different sets of quenched data with /3 = 6.0, 
6.2 and 6.4 and using either the Wilson or the 
"improved" SW-Clover action. The parameters 
used in each simulation are listed in Table 1 and 
the main results we have obtained are reported in 
Tables ^ and ||. The data at /? = 6.4 have been 
used only for an exploratory study. The physical 
volume and the time extension of the lattice are 
too small to be considered reliable. 
The MS mass values in Table || show that the 
non-perturbative approach allows a determina- 
tion of the quark masses with a smaller error 
than the standard perturbative approach. In 
the /3 range studied, there is no statistical evi- 
dence for an " a" dependence of the quark masses 
for the Clover action. For the Wilson action, 
a mild tendency in the quark masses to de- 
crease with increasing [3 exists but an extrapo- 
lation to the continuum limit is not reliable. The 
MS mass values in Tables || and ^ show a very 
good agreement between the quark masses ex- 
tracted from the spectroscopy and from the Ward 



Identity using the non perturbative determina- 
tions of the renormalization constants, while the 
same comparison using the perturbative values 
of Za and Zp is much poorer We be- 

lieve that the best estimates for the strange and 
charm quark masses are obtained from the non- 
perturbative approach using Wilson and Clover 
data aX (3 — 6.0 and (3 = 6.2. Our best results 
are m^{2 GeV) = (123 ± 4 ± 15) MeV and 
m|^(2 GeV) = (1525 ± 40 ± 100) MeV. The 
first error is the statistical one and the second 
is the systematic error estimated extracting the 
quark masses from different mesons The val- 
ues of quark masses we have obtained are in good 
agreement with previous determinations |p],p|-p^ . 
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